Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN I*TeX style file vl.4) 

Brown Dwarfs and the Cataclysmic Variable Period 
Minimum 

U. Kolb^^ and I. Baraffe^ 

^ Astronomy Group, University of Leicester, Leicester, LEI 7RH 

2 Ecole Normale Superieure de Lyon, C.R.A.L. (UMR 5574 CNRS), F-69364 Lyon Cedex 07, France 



1—5 

oo 
(N 



> 

oo 

o 
(^ 

6 



June 1999 



% 



ABSTRACT 

Using improved, up-to-date stellar input physics tested against observations of low- 
mass stars and brown dwarfs we calculate the secular evolution of low-mass donor cat- 
aclysmic variables (CVs), including those which form with a brown dwarf donor. Our 
models confirm the mismatch between the calculated minimum period (Pmin — 70 min) 
and the observed short-period cut-off (~ 80 min) in the CV period histogram. We 
find that tidal and rotational corrections applied to the one-dimensional stellar struc- 
ture equations have no significant effect on the period minimum. Theoretical period 
distributions synthesized from our model sequences always show an accumulation of 
systems at the minimum period, a feature absent from the observed distribution. We 
suggest that non-magnetic CVs become unobservable as they are effectively trapped 
in permanent quiescence before they reach Pmin, and that small-number statistics may 
hide the period spike for magnetic CVs. 



Key words: accretion, accretion discs — binaries: close - 
ables — stars: evolution — stars: low-mass, brown dwarfs. 
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1 INTRODUCTION 

, Cataclysmic variables (CVs) are semi-detached binaries 

' with a white dwarf (WD) primary and a low-mass main- 
sequence companion that transfers mass to the WD through 

, Roche-lobe overflow (e.g. Warner 1995) . Orbital periods are 
known for about 400 systems. Two marked features stand 

\ out in the orbital period distribution: a dearth of systems 
in the range 2-3 hr, usually referred to as the "period gap", 
and a sharp short-period cut-off at ~ 80 min, the "period 
minimum" Pmin (e.g. Ritter & Kolb 1998). J0132-6554 at 
P = 77.8 min currently marks the short end of the bulk 
of the distribution, while the single system V485 Cen at 
P — 59 min is the notable exception. Six AM CVn-type 
CVs with yet shorter periods are interpreted as CVs with 
helium donors. 

In a semi-detached binary the Roche lobe filling star's 
mean density p determines the orbital period P ahnost 
uniquely (e.g. King 1988), Ph = k/pQ (where k ~ 8.85 is 
only a weak function of the mass ratio, P^ is the period in hr, 
Pq the mean density in solar units) . The density increases for 
donors evolving along the hydrogen-burning main sequence 
towards smaller mass. Approaching the hydrogen-burning 
minimum mass the increasing electron degeneracy induces 
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structural changes such that further mass loss reduces p. 
Hence during the donor's transition from a main-sequence 
star to a brown dwarf (BD) the secular mean orbital period 
derivative changes from negative to positive. This "period 
bounce" has long been identified with the minimum period 
of CVs (Paczyiiski 1981, Paczyhski & Sienkiewicz 1981; Rap- 
paport. Joss & Webbink 1982; D'Antona & Mazzitelh 1982). 

The actual period Ptum at which a CV bounces depends 
on the ratio r = t-Ku/tM of the secondary's thermal time 
fxH ~ GM^/RL and mass transfer time Im = M/(— M2), 
the characteristic timescales restoring and perturbing equi- 
librium. If r is small Pmin is short. Hence Ptum is sensitive 
both to the orbital angular momentum loss rate (which de- 
termines the transfer rate —M2) and to the interior structure 
of the donor. Paczynski (1981) was the first to point out that 
Ptum is close to 80 min if gravitational wave radiation drives 
the mass transfer. Since then, models with different input 
physics have been employed to verify a quantitative agree- 
ment between the observed Pmin and the calculated Ptum- 
Recent calculations (Kolb & Ritter 1992; Howell, Rappa- 
port & Politano 1997) notoriously give Ptum too short by 
typically 10%. 

Most of these evolutionary calculations are based on in- 
put physics which is rather approximate for very low mass 
(VLM) objects, i.e. donors with mass ;$ IM©. Because of 
their relatively high central densities and low central temper- 
atures, correlation effects between particles dominate and 
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must be taken into account in the equation of state (cf. 
Chabrier & Baraffe 1997, and references therein). Below 
about Teff ^ 4000 K (M <, O.6M0), molecules become stable 
and dominate the atmospheric opacity, being responsible for 
strong non-grey effects and significant departure of the spec- 
tral energy distribution from blackbody emission (cf. AUard 
et al. 1997, and references therein). 

Because of significant efforts devoted to the complex 
physics of low-mass stellar/substellar objects, the theory 
describing them has improved considerably in the past few 
years. As a result, the latest generation of stellar models for 
low-mass stars and brown dwarfs is now able to reproduce 
observed properties of field M-dwarfs with unprecedented 
accuracy (cf. AUard et al. 1996; Marley et al. 1996; Chabrier 
& Baraffe 1997; Burrows et al. 1997; Baraffe et al. 1998). 

In this paper we use the Baraffe et al. models (1995, 
1997, 1998, henceforth summarized as BCAH) — briefly re- 
viewed in Sect. H — to calculate the secular evolution of CVs 
in the vicinity of Ptum (Sect. 3.1). We test if tidally distorted 



stellar models lead to a significant increase of Ptum over the 
value for spherical stars, as claimed by Nelson et al. 1985 
(Sect. 3.2), and address the problem of the missing "period 
spike" at Pmin (Sect. 3.3). This predicted accumulation of 
systems at Ptum is caused by the slow velocity in period 
space close to period bounce which increases the detection 
probability. The spike has been a dominant feature in syn- 
thesized period distributions obtained from theoretical CV 
population models (Kolb 1993; Kolb et al. 1998), yet is ab- 
sent in the observed CV period distribution. In these earlier 
population models the donor star was approximated as a 
(bi-)polytrope, and CVs forming with a donor mass smaller 
than 0.10 M© were not considered. Here we investigate if 
the period spike persists in period distributions synthesized 
from full BCAH models. For the first time, we include CVs 
which form already with BD donors, i.e. systems which did 
not evolve through period bounce. In Sect, q alternative ex- 
planations for the missing period spike are discussed. In par- 
ticular, we focus on the difference between magnetic and 
non-magnetic CVs. 



2 THE SECONDARY'S INTERNAL 
STRUCTURE 

Baraffe et al. (1998, and references therein) give a brief ac- 
count of the input physics used in the most recent low-mass 
star models we apply here. The main strengths of these mod- 
els are in two areas: the microphysics determining the equa- 
tion of state (EOS) in the stellar interior, and the outer 
boundary condition based on non-grey atmosphere models. 
The models employ the Saumon, Chabrier and Van 
Horn (1995) EOS which is specifically calculated for VLM 
stars, BDs and giant planets. The EOS is an important in- 
gredient for our analysis, since it determines the mechanical 
structure of the donor stars, and thus their mass-radius rela- 
tion. Models based on this EOS have been successfully tested 
against stars in detached eclipsing binary systems (Chabrier 
& Baraffe 1995) and field M-dwarfs (cf. Beuermann et al. 
1998). The Saumon et al. (1995) EOS has been successfully 
compared with recent laser-driven shock wave experiments 
performed at Livermore. These probe the complex regime of 
pressure dissociation and ionization which is so characteris- 



tic for the objects we are studying here (cf. Saumon et al. 
1998). 

The second important ingredient of our VLM star and 
BD models concerns the outer boundary condition, which 
determines the thermal properties of an object and thus the 
mass-effective temperature (Tcs) relation. As demonstrated 
by Chabrier & Baraffe (1997), evolutionary models with a 
grey atmosphere outer boundary, e.g. the standard Edding- 
ton approximation, overestimate Tcb for a given mass, and 
yield too large a minimum hydrogen burning mass (MHBM) . 
Because of tremendous recent progress in the field of cool at- 
mosphere models (see e.g. the review of AUard et al. 1997) 
these now provide realistic atmosphere profiles and synthetic 
spectra which we use as a more realistic outer boundary con- 
dition. 

Several observational tests confirm the success of 
evolutionary models based on these improvements, e.g. 
mass-magnitude relationships, colour-magnitude diagrams 
(Baraffe et al. 1997, 1998), mass-spectral type relationships 
(Baraffe & Chabrier 1996), and the first cool BD GL 229B 
(AUard et al. 1996). Although some discrepancies between 
models and observations remain (see Baraffe et al. 1998), 
uncertainties attributed to the input physics now seem sig- 
nificantly reduced. 

We employ these one-dimensional BCAH models in the 
context of CVs within the usual Roche approximation. The 
mass transfer rate is calculated as a function of the dif- 
ference between donor radius and Roche radius, following 
Ritter (1988). In a separate set of calculations we apply 
tidal and rotational corrections to the one-dimension stellar 
structure equations (Chan & Chan 1979) to account for the 
non-spherical shape of Roche equipotential surfaces. This 
allows us to test the results by Nelson et al. (1985), who 
found that their distorted stellar models give a value for 
Ptum which is about ~ 10% longer than for their spherical 
stellar models. 



3 MODEL CALCULATIONS 

3.1 Evolutionary sequences 

Using the BCAH stellar evolution code we recalculated the 
secular evolution of short-period CVs for various initial com- 
ponent masses. We assume that orbital angular momentum 
is lost by gravitational wave emission (e.g. Landau & Lif- 
shitz 1958) and via an isotropic stellar wind from the WD. 
The wind removes the accreted mass with the WD's specific 
orbital angular momentum from the binary. Hence the WD 
mass is constant throughout the evolution. 

As we focus on the period minimum we do not consider 
the evolution at periods longer than 2-3 hr. In a second 
paper (Baraffe & Kolb 1999; see also Kolb & Baraffe 1999) 
we discuss properties of such sequences with more massive 
donors. They depend on the rather uncertain strength of the 
magnetic braking thought to dominate the evolution of CVs 
above the period gap (see e.g. King 1988 for a review). 

We calculated three sets of evolutionary sequences. In 
Set A the WD mass is 0.6 M©, while the initial donor masses 
range from 0.27 M© to 0.04 M© (see Fig. |). Set B is similar 
but for a WD mass 1.0 Mq. In Set C the initial donor mass 
is fixed at 0.21 Mq , while the WD masses range from 1.2 Mq 
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Table 1. Model parameters of sequences in Set C (initial donor 
mass 0.21 Mq; age of donor at birth of CV 0.61 Gyr). Mi is the 
WD mass, tf the age since formation of the CV when Tgff = 
900 K, Mturn is the donor mass at period bounce, mturn the 
transfer rate at period bounce (in 10~-'^^Mq yr~^). 

Mi{Mq) tf (Gyr) Ptum (hr) Mtuvn (Mq) rhtuvn 



0.30 
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0.40 
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0.70 
0.80 
1.00 
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6.35 
5.94 
5.61 
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Table 2. Characteristic quantities along the sequence with initial 
donor mass 0.21 Mq and WD mass 0.6 Mq. t is the age in Gyr, 
P^ the orbital period in hours, M the donor mass in Mq ; T^fj the 
effective temperature in K; L corresponds to log L/Lq; R is the 
radius in solar units, and M the transfer rate in MQyr"-"^. 



t (Gyrs) 


Ph 


M 


^off 


L 


R/Rq 


logM 


0.00 


2.080 


0.2103 


3311. 


-2.272 


0.2242 


-16.00 


0.25 


2.038 


0.1975 


3288. 


-2.313 


0.2168 


-10.34 


0.50 


1.970 


0.1862 


3264. 


-2.363 


0.2077 


-10.35 


0.75 


1.902 


0.1753 


3239. 


-2.415 


0.1987 


-10.37 


1.00 


1.832 


0.1648 


3212. 


-2.470 


0.1898 


-10.38 


1.25 


1.764 


0.1545 


3177. 


-2.530 


0.1810 


-10.39 


1.50 


1.696 


0.1446 


3133. 


-2.596 


0.1725 


-10.41 


1.75 


1.628 


0.1351 


3087. 


-2.665 


0.1640 


-10.43 


2.00 


1.556 


0.1259 


3044. 


-2.736 


0.1554 


-10.44 


2.25 


1.483 


0.1169 


2993. 


-2.815 


0.1469 


-10.45 


2.50 


1.409 


0.1082 


2929. 


-2.904 


0.1384 


-10.46 


2.75 


1.336 


0.0996 


2849. 


-3.007 


0.1300 


-10.47 


3.00 


1.265 


0.0912 


2743. 


-3.129 


0.1217 


-10.48 


3.25 


1.200 


0.0830 


2607. 


-3.275 


0.1140 


-10.49 


3.50 


1.148 


0.0750 


2423. 


-3.456 


0.1071 


-10.50 


3.75 


1.119 


0.0675 


2205. 


-3.665 


0.1017 


-10.54 


4.00 


1.114 


0.0607 


1972. 


-3.890 


0.0980 


-10.60 


4.25 


1.130 


0.0549 


1750. 


-4.117 


0.0958 


-10.68 


4.50 


1.157 


0.0501 


1570. 


-4.317 


0.0946 


-10.77 


4.75 


1.189 


0.0463 


1417. 


-4.503 


0.0939 


-10.86 


5.00 


1.223 


0.0432 


1280. 


-4.682 


0.0935 


-10.94 


5.25 


1.257 


0.0406 


1176. 


-4.830 


0.0934 


-11.02 


5.50 


1.290 


0.0384 


1091. 


-4.961 


0.0934 


-11.10 


5.75 


1.321 


0.0365 


1025. 


-5.070 


0.0934 


-11.17 


6.00 


1.350 


0.0350 


968. 


-5.168 


0.0934 


-11.23 


6.25 


1.377 


0.0336 


918. 


-5.260 


0.0934 


-11.29 


6.35 


1.387 


0.0331 


900. 


-5.295 


0.0935 


-11.31 



Table 3. Same as Table H, but for a BD CV sequence with initial 
donor mass M2 = 0.07 Mq and WD mass 0.6 Mq. 



Figure 1. Mass transfer rate versus orbital period for sequences 
of Set A (white dwarf mass 0.6 Mq). Labels indicate the initial 
donor mass (in Mq). 



to 0.3 Mq; see Tab. hi. At turn-on of mass transfer the sec- 
ondary was either on the ZAMS {Ah > 0.085 Mq), or had 
an age of 2 Gyrs (M2 < 0.085 Mq). All donors are hydrogen- 
rich and have solar metallicity {X = 0.70, Z = 0.02). The 
sequences were terminated when the donor's effective tem- 
perature TcH was smaller than 900 K. This is the lower limit 
of the temperature range covered by the present non-grey at- 
mosphere models (Hauschildt et al. 1999) available for this 
study. 

Some quantities of interest for a typical sequence start- 
ing immediately below the periodgap are shown in Fig. Q 
and given in Tab. H (see also Tab. fa for a sequence initiating 
mass transfer from a BD donor). 

Fig. n] confirms the well known effect that systems with 
different initial donor masses rather quickly join a uniform 
evolutionary track (Stehle et al. 1996). Most systems of 
Set A undergo period bounce at Ptum — 67 min, which is 
only slightly longer than the corresponding Ptum ~ 65 min 



t (Gyrs) Ph 



M 



R/Rq log M 



0.00 
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2.00 


1.358 


0.0347 


943. 


-5.214 


0.0935 


-11.24 


2.01 


1.359 


0.0346 


940. 


-5.219 


0.0935 


-11.24 



found with Mazzitelli's models (Kolb & Ritter 1992; see also 
the dotted curve in Fig 0, lower panel). Note that a grey 
outer boundary condition like the standard Eddington ap- 
proximation yields a smaller radius for given donor mass 
M2 < 0.2Mq, and a higher MHBM. A test calculation with 
fixed mass loss rate and the Eddington approximation gives 
Pturn shorter by 5 min than for the corresponding sequence 
based on non-grey atmosphere models. Note also that CVs 
forming with fairly old and massive brown dwarf donors (age 
<; 2 Gyr, mass 0.05 — 0.07Mq) would populate the period 
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Figure 2. Time evolution of the sequence with O.6M0 WD mass, 
0.21Mq initial donor mass sequence. The time t at turn-on of 
mass transfer has been set arbitrarily to 1 Gyr. Upper panel: 
Effective temperature Tj.g (dash-dotted; left-hand scale) and ra- 
dius r2 (dashed), mass m2 (solid), both in solar units, right— hand 
scale, of the donor. Lower panel: Mass transfer rate rh (in 10"^*^ 
Mq yr~^; solid) and orbital period P (dashed). 



regime shortwards of Ptum. V485 Cen might be such a sys- 
tem. 

For Set B we have Ptum ~ 69 min., slightly longer than 
for Set A since gravitational radiation losses are higher for 
larger WD mass, hence the mass transfer rate is larger. The 
dependence of Ptum on WD mass is most easily seen in 
Fig. 0, upper panel. 

The effective temperature of the donor as a function of 
period is shown in Fig. a The asymptotic convergence of 
sequences with different initial donor mass and the insensi- 
tivity of M to the WD mass result in an essentially unique 
relation TeS (P) on the non-degenerate branch of the track 
(before period bounce occurrs). At Ptum, and on the de- 
generate branch (after period bounce) there is a significant 
spread in Tcff for given P. In the same diagram dots indicate 
the claimed location of AL Com (Howell et al. 1998), WZ 
Sge, TY Psc, V592 Cas and HU Aqr (Ciardi et al. 1998). 
For these systems the effective temperature of the donor 
has been derived via spectral fittings in the near-infrared 
using the same atmosphere models (Hauschildt et al. 1999) 
that serve as outer boundary condition of our stellar mod- 
els. Since the near-IR emission of CVs is produced by both 
the secondary and the accretion disc, the estimate of the 
secondary's contribution may be rather uncertain because 



Figure 3. Effective temperature versus orbital period for vari- 
ous evolutionary sequences. Upper panel: Sequences with different 
WD masses. Solid: 1.2M0, O.6M0 O.SM©, in order of decreas- 
ing Ptum. Dashed: O.7OM0, but with a low-nictallicity donor 
(Z = 0.006). Dots show the claimed location of AL Com (Howell 
et al. 1998), WZ Sge, TY Psc, V592 Cas, and HU Aqr (Ciardi et 
al. 1998). Lower panel: Sequences (of Set A) with different ini- 
tial donor masses (solid; see also Fig. hi). Dotted: a 0.7Mq WD 
mass sequence, calculated with Mazzitelli's (1989) stellar evolu- 
tion code. 



of the lack of reliable accretion disc models. We thus do 
not regard the discrepancy between model calculations and 
observational data points (cf. Fig.H) as significant. 



3.2 Sequences with non— spherical donors 

We implemented the rotational correction scheme of Chan 
& Chau (1979) to investigate the effect of rotation and tidal 
distortion on the donor star. Assuming solid-body rotation, 
we recalculated a sequence with O.2IM0 initial donor mass, 
I.OM0 WD mass and J = Jgr. We find that Ptum is hardly 
affected; it is longer by only 1 min compared to the cor- 
responding sequence with spherical models (from Set B). 
We have performed several other experiments, with differ- 
ent initial secondary masses and mass transfer rates (con- 
stant, or as given by gravitational radiation), including the 
case considered by Nelson et al. (1985), i.e. Mi = I.OM0, 
M2 = O.4M0 , J = Jgr. In none of these sequences we could 
reproduce the significant effect (~ 10%) rotational and tidal 
corrections had on Ptum in the calculations by Nelson et al. 
(1985) based on the same scheme of Chan and Chau (1979). 
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We do not think that the discrepancy between Nelson 
et al. (1985) and our work results from differences in the 
EOS and opacities. Indeed, although our input physics dif- 
fers from those adopted by Chan & Chan (1979), we find 
the same quantitative differences quoted by these authors 
between the properties {L, R and Tets) of spherical and ro- 
tating ZAMS stars (cf. their Table 1). We also reproduce 
the effects quoted in the very recent paper by Sills & Pin- 
sonneault (1999). They considered a 0.7 Mq ZAMS star 
with equatorial rotational velocity 145 km s^^ and find a 
slightly increased surface luminosity and effective temper- 
ature (AlogL = 3.2%, ATeff = 100 K) compared to non- 
rotating models. For the same case we find AlogL = 4.5%, 
ATeff = 90 K. 

We thus are confident in our results and conclude that 
rotational and tidal corrections as described by Chan and 
Chan (1979) can be neglected for CVs, even when the sys- 
tems are close to period bounce. Note that, although the 
secondary spins up in evolving toward Ptum, the key quan- 
tity defined by Chan and Chau (1979), 

2a;2i?3 
a = -- 



3 GMr ' 

which appears in the momentum equation and represents 
the ratio of rotational to gravitational potential gradient, 
remains essentially constant at the donor star surface during 
the evolution (recall that the angular velocity lo'^ oc P~^ cx 
M2/i?2)- The maximum value of a, reached at the surface, 
never exceeds 0.07. 

Given the small differences between spherical and 
tidally distorted models we use the less cpu-time intensive 
calculations with spherical stars to derive a period distribu- 
tion and analyse its properties close to Pmin. 



3.3 Generating a period histogram 



When the evolutionary sequences presented in Sect. 3.1 are 



convolved with an appropriate CV formation rate they give 
a theoretically predicted period histogram in the vicinity 
of Pmin (see Kolb 1993 for a detailed description of the 
population synthesis technique). De Kool (1992) and Poli- 
tano (1996) calculated the formation of CVs with standard 
assumptions about common envelope evolution and mag- 
netic braking (see e.g. Kolb 1996 for a review). We com- 
puted population models for all these formation rate mod- 
els; in addition, to show the differential effect of the main 
parameters more clearly, we considered subpopulations for 
a given WD mass with a time-independent formation rate 
6(logM2) = const, in certain logM2 intervals. 

In Figs. 0, we show period distributions for two differ- 
ent CV subpopulations — PI and P2 — which allow one to 
study the main effects that determine the structure of the 
histogram at Pmin. The histograms are for either a volume- 
limited sample, or a sample where individual systems have 
been weighted by M" , for various values of a. This mimics 
selection effects which affect the observed sample (cf. Kolb 
1996). The subpopulations formally correspond to a Galac- 
tic disc age of 6 Gyr. In a somewhat older population the 
edge of the volume-limited distributions at ~ 1.4 hr would 
appear at slightly longer P, while the M" -weighted cases 
would hardly change. 



PI. Period histograms obtained from Set A (Fig- 
ure M). These correspond to the period distribution of a sub- 
set of CVs with WD mass O.6M0 (assuming that this mass 
does not change through the evolution). The same analysis 
has been performed for sequences of Set B, with very similar 
results, which we therefore do not show. 

As noted above we do not consider the evolution of CVs 
with donors > O.2IM0, i.e. of CVs that form above the pe- 
riod gap. In the disrupted magnetic braking model for the 
period gap (e.g. King 1988) these systems would all reap- 
pear ^ 10* yr after formation at the lower edge of the gap, 
with a donor in thermal equilibrium. A BCAH stellar model 
in thermal equilibrium with mass ~ O.2IM0 would fill its 
Roche lobe at P = 2.1 hr, the observed lower edge of the 
period gap. Hence we simulate the contribution of systems 
that form above the gap to the period distribution below 
the gap by increasing the formation rate in a narrow bin in 
M2-space at 0.21 Mq by a large factor. 

For the first time, we explicitly include systems which 
form with a brown dwarf donor star. The formation rate of 
such BD CVs is not known. Survival of the common enve- 
lope phase is crucial and could be a problem as the maxi- 
mum orbital energy available to eject the envelope, roughly 
oc M2, is small. Simulations by Politano (1998, priv. comm.) 
show that the formation of BD CVs is possible when the 
same formalism is applied as for more massive donor stars. 
For our purposes we just extrapolate the birthrate function 
&2(logM2) down to the smallest initial donor mass (O.O4M0) 
considered here. 

In particular, we use 



&2(logM2) 




0.207 < Afa < 0.210 
0.090 < Ma/Me < 0.207 
0.040 < Mz/Mq < 0.090 



(1) 



with various combinations of the free parameters fci and fca, 
describing the relative weight of systems forming above the 
period gap and with a brown dwarf donor, respectively. The 
numerical factor in front of fci is chosen such that for fci = 1 
as many CVs form with donor masses 0.21 < M2/M0 < 1.0 
("above the gap") as with masses 0.09 < M2/M0 < 0.21 
( "below the gap" ) . The results are not sensitive to the choice 
of the boundaries at O.2O7M0 and O.O9M0. So, chiefly, fci is 
the ratio of systems forming above the gap to systems form- 
ing below the gap with a non-degenerate donor, while ^2 is 
the ratio of the formation rate of CVs with degenerate and 
non-degenerate donors. Figure q shows period histograms 
for various combinations of fci, ^2 and a. 

P2. A period histogram obtained from Set C, with 
the full time-dependent formation rate calculated by de 
Kool (1992; his model 3). This distribution, shown in the 
top panel of Fig. 0, illustrates the fine structure of the pe- 
riod spike at Ptum as a result of the different WD masses in 
the sample. To obtain this histogram, we approximated se- 
quences with an arbitrary initial donor mass M2 < 0.21 M0 
by the 0.21 M0 sequence from M2 onwards. This introduces 
a slight error at the onset of mass transfer and for sequences 
which would form close to Ptum. An explicit comparison 
of results from this simplified method for O.6M0 WD mass 
with the distribution obtained from the full set of sequences 
(set A) shows that these deviations are negligible as long 
as ^2 ^ 2. To estimate the contribution from systems that 
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Figure 4. Orbital period distributions n(logP) for CVs with a 0.6 Mq WD, synthesized from Set A and formation rate (hh, arbitrarily 
normalised, for different fractions fci of systems forming above the period gap, fractions ki of CVs forming with BD donors, and 
detectability parameters a (see text for details). The solid curves are the corresponding cumulative distributions N(> P). Left: k2 = 1, 
a = 1 in all models, fci as labelled. Right: fci = 1 in all models, k2 and a as labelled. 



form above the period gap the actually calculated evolution- 
ary tracks have been extended to donor masses > O.2IM0 
by assuming a simple main-sequence mass-radius relation 
(such that P oc A/2) and a constant mass transfer rate 
(2 X 10-« Moyr-i). 

The theoretical period distributions in Figs. ^ and q show 
that the collective period minimum Pmin coincides with the 
bounce period Ptum of individual evolutionary sequences, 
and is about 10 min. shorter than the observed period min- 
imum (e.g. Ritter & Kolb 1998). This confirms tentative 
conclusions reached earlier by Kolb (1993) who considered 
populations constructed from simplified polytropic stellar 
models with artificial outer boundary conditions, and re- 
stricted to initial donor masses > O.O9M0. 

The subpopulations for WD masses 0.6 M© and 1.0 Mq 
show that the period spike does not disappear or broaden 
significantly when CVs with brown dwarf donor stars are 
included in the population. This is true unless the majority 
of newly forming CVs actually are BD CVs (^2 ^ 5). But if 
BD CVs were dominant the bulk of systems would populate 
the period regime shortwards of Ptum, making this a fairly 
unlikely possibility in view of the observed value of Pmin. 

In fact, in population models which do not explicitly 



overemphasize BD CVs (^2 ^ 2) the effect of these sequences 
on the overall shape of the period spike is generally small. 
A fair representation of the period spike is obtained by the 
simple procedure used to generate the period distribution 
P2. The double spike at Pmin in the P2 histogram (Fig. pi 
upper panel) is a result of the double-humped WD mass 
distribution (see e.g. de Kool 1992): high-mass C/O WD 
systems cluster in the longer-period spike, low-mass He- 
WD systems in the shorter-period spike. This fine structure 
of course disappears altogether in small samples. 

The overall shape of the weighted distribution is not 
sensitive to moderate a; there is hardly any difference be- 
tween histograms calculated with a between 0.5 and 1.5. 
Increasing a further tends to decrease the amplitude of the 
spike. We obtain a truly flat distribution, similar to the ob- 
served one (Fig. H), only for q <; 6. 



4 DISCUSSION 

The results presented in the previous section show that 
despite the signiflcant improvement of low-mass star and 
brown-dwarf models, the period histogram for short-period 
CVs synthesized from these models still appears to be in 
conflict with observations. If the evolution of short-period 
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Figure 5. Upper panel: synthesized orbital period distribution 
P2 (see text), for a = 1.5. Middle and lower panel: two arbitrary 
samples of 33 systems drawn from the distribution in the upper 
panel. 



Figure 6. The observed orbital period distribution of dwarf novae 
(upper panel) and polars (lower panel), for systems with P^ < 2.1. 
The solid curve is the normalized cumulative distribution. Data 
from Rittcr & Kolb (1998). 



CVs is driven by the emission of gravitational waves, period 
bounce occurs at ~ 70 min for most systems, resulting in 
a collective minimum period Pmin which is 10 min shorter 
than the observed one. The contribution of CVs forming 
with BD donors is either negligible or would cause Pmin to 
be yet shorter. As a result of period bounce, the simulated 
period histograms show an accumulation of systems near 
Pmin, the period spike. The only way to remove the period 
spike effectively from the distribution is to impose a very 
steep dependence of the system's detectability d on one of 
the evolutionary parameters, e.g. oc M® or steeper. 

The catalogue of Ritter & Kolb (1998) hsts 129 CVs 
with periods < 2.10 hr. Among these are 33 polars, 78 dwarf 
novae (excluding intermediate polars), and 6 AM CVn sys- 
tems. The remaining 12 systems are either intermediate po- 
lars, novae, novalikes, or of unknown type. Hence essentially 
all known short-period CVs are either dwarf novae, i.e. sys- 
tems with an unstable accretion disc, or discless polars. 



4.1 The missing period spike 

Almost all dwarf novae are detected in outburst. It is tempt- 
ing to identify the accretion disc and its outburst properties 
as the cause for a strong M-dependent selection effect. A few 
short-period dwarf novae, sometimes referred to as WZ Sge 
stars, have very long outburst recurrence times irec, the most 
extreme example being WZ Sge itself with tree — 30 yr. It 



has long been noted that low-M CVs might have escaped de- 
tection if their outburst interval is significantly longer than 
the period since beginning of systematic monitoring and sur- 
veying of the sky with modern means - a few decades. For 
long iree the relative detection probability of dwarf novae 
scales as d oc 1/tree, suggesting that iree oc M^® or steeper. 
In practice this very steep dependence would require that 
tree ^ oo for M ^ few X lO""M0yr"\ i.e. low-M CVs 
would not undergo outbursts at all. 

Various plausible physical models which naturally ac- 
count for this property have been suggested. Evaporation 
of accreted material into a hot corona could prevent the 
disc from accumulating the critical surface mass density re- 
quired to launch a heating wave (Meyer-Hofmeister et al. 
1998; see also Liu et al. 1997). Systems in such a permanent 
quiescence would be optically very faint but should still emit 
about 10% of their accretion luminosity Lace = GMiM / R\ 
{Mi,Ri is the WD's mass and radius) in X-rays. Alterna- 
tively, King (1999) argues that systems approaching the pe- 
riod minimum become magnetic ejectors even if the WD's 
magnetic field is weak. In this case an accretion disc can 
no longer form, and the systems become unobservable. WZ 
Sge, where the WD is suspected to be weakly magnetic (see 
e.g. Lasota et al. 1999, and references therein), would be a 
marginal disc accretor (Wynn & Leach 1999). 

Although this can explain the absence of a period spike 
at Pmin for dwarf novae, it certainly does not apply to the 
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discless polars. An independent observational selection ef- 
fect operating on polars in the same period/mass transfer 
rate range and with the same net result as a steep increase 
of tree for dwarf novae seems a highly unlikely coincidence 
(but see Meyer & Meyer-Hofmeister 1999). In other words: 
polars should show a period spike, even though dwarf novae 
do not. Could small-number statistics hide the period spike 
for polars? To investigate this we performed a number of 
Monte Carlo experiments, where we drew a sample of either 
33 or 78 systems from an underlying period distribution like 
the one shown in Fig. |5| (upper panel). As the middle and 
lower panels of this figure show, the smaller samples give 
a surprisingly wide variety of distributions, with typically 
20% — 25% of them showing no sign of a period spike at all. 
In contrast, in the larger sample the spike is prominent in 
more than 90% of all cases. A KS test confirms this impres- 
sion: We wish to quantify the difference between model and 
observed distribution that is due to a different morphological 
shape. Therefore we have to exclude effects arising from dif- 
ferent values for the lower edge Pgap of the period gap and 
for Pmin — the latter effect is considered in Sect. 4.2. To 



achieve this we rescale the calculated distribution such that 
it matches the range of the observed distribution before per- 
forming the KS test. Specifically, the period axis C = logP 
of the theoretical distribution is rescaled according to 






(2) 



where Ci — log Pmin, Cu = log Pgap denote the calculated 
period minimum and lower edge of the gap, and Oi, Ou 
the corresponding observed values. We applied (g) to the 
model distribution shown in Fig. g (Pmin ~ 64 min, Pgap = 
124 min), assuming the observed values Pmin = 78 min, 
Pgap = 113 min. The maximum significance level for reject- 
ing the null hypothesis that the observed sample is drawn 
from this rescaled model distribution is 0.34 for polars, 0.98 
for dwarf novae. If we use Pgap = 130 min for the observed 
lower egde of the gap, the rejection significance for polars 
rises to 0.91. This is solely due to the so-called "114 min 
spike", see e.g. Ritter & Kolb 1992, still a significant feature 
in the observed distribution. 



4.2 The mismatch between observed and 
calculated Pmin 

Neither the small number statistics for polars nor the sug- 
gested detectability function for dwarf nova can make the 
minimum period of the observed distribution significantly 
longer than in the underlying intrinsic distribution. To 
achieve this by the latter effect, dwarf novae have to become 
unobservable long before they reach Pturn. With M as the 
most likely control parameter determining the detectability 
this is difficult to achieve, as M is almost constant on the 
non-degenerate branch above Pturn. (This problem would 
be less severe if the main control parameter were the donor 
mass, the mass ratio, or the orbital separation. A particu- 
larly steep dependence of the discovery probability on any 
of these could be achieved if mass transfer cycles, similar to 
the ones discussed by King et al. 1996, 1997, existed in CVs 
below the period gap, but changed their character discon- 
tinuously before Pturn.) 

To test the effect of small number statistics we drew 



1000 samples of either 33 or 78 systems from the theoret- 
ical distribution shown in Fig. and registered the short- 
est period Ps of each sample. The parent distribution has 
a short-period cut-off at 64.4 min. We found that 99.9% 
of the larger samples have Ps < 67.0 min (99.0% have 
Ps < 66.6 min, 90.0% have P^ < 65.6 min). The median 
is Ps = 65 min, the longest Ps we found is 67.03 min. Sim- 
ilarly, 99.9% of the smaller samples have Ps < 71.0 min 
(99.0% have Ps < 68.2 min, 90.0% have Ps < 66.8 min), 
with a median Ps = 65.4 min, and 71.2 min as the longest 
Ps- The rare cases where Ps was as long as 70 min still fall 
well short of the observed value Pmin = 78 min. 

Given this, it might be that the Pmin mismatch is caused 
by evolutionary effects after all — effects not accounted for 
in our models. 

The chemical composition of the secondary affects the 
value of Pturn, by changing the parameter of the donor at a 
given mass, hence changing both Im (via the angular mo- 
mentum loss time) and tKH. Stehle et al. (1997) pointed out 
that Pturn is slightly shorter for CVs with low-metallicity 
donors. They found dPtum/dlogZ = 0.084 hr. Our calcula- 
tions confirm this (Ptum = 67.41 mins. for Z = 0.02, while 
Pt^,„ = 64.70 mins for Z = 0.006; both with Mi = 0.70Mq). 
Although Pturn increases for higher metallicity donors the ef- 
fect is much too small to account for a mismatch of 10 mins. 
It has been noted that a larger than expected fraction of CVs 
above the gap have a nuclear-evolved secondary (Beuer- 
mann et al. 1998, Baraflte & Kolb 1999, Kolb & Baraflte 
1999). When these secondaries become fully convective their 
helium content is ~ 0.5, giving a Ptum which is in fact 
~ 7 mins. shorter than for hydrogen-rich donors. 

Residual shortcomings in the EOS and the atmosphere 
profiles cannot yet be ruled out as the cause for the Pmin 
discrepancy. A quantitative estimate of uncertainties in the 
mass-radius relation from the treatment of the EOS (see 
Saumon et al. 1995) is difficult. The most profound observa- 
tional test is against stellar parameters determined for com- 
ponents in (detached) eclipsing binaries; yet to date there 
are no such systems with M2 < O.2M0. The region near 
Pturn involves effective temperatures lower than 2600 K (cf. 
Fig. pi), below which grains form. These affect both the at- 
mosphere spectrum and profile. Although the atmosphere 
models used for this study are grainless, calculations based 
on preliminary atmosphere models including the formation 
and absorption of grains (Allard 1999; Baraffe & Chabrier 
1999) do not yield a longer Ptum. We thus do not expect an 
improvement of the situation with the forthcoming genera- 
tion of dusty atmosphere models. 

The Pmin mismatch could also be due to uncertainties 
in the calculated value of Ptum inherent to the very concept 
of the Roche model. Strictly valid only for point masses, 
its applicability to extended donors relies on the fact that 
the stars are usually sufficiently centrally condensed. This 
is not necessarily a good approximation for fully convec- 
tive stars which are essentially polytropes of index n — 3/2. 
We note that Uryu & Eriguchi (1999) considered stationary 
configurations of fiuids describing a synchronously rotating 
polytropic star in a binary with a point mass companion. For 
n — 3/2 polytropes they find a Roche radius that is typically 
1 — 2% smaller than estimated from Eggleton's (1983) ap- 
proximation for a given orbital separation, but at the same 
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time ~ 4% larger than the radius of a non-rotating poly- 
trope with the same mass. 

Alternatively, an obvious way to increase Pmin is to in- 
crease the orbital angular momentum loss rate J over the 
value Jgr set by gravitational radiation. We find Pmin — 
83 min (up from 69 min) for J = 4 x Jgr and Mi = IM0, 
a much smaller increase than quoted by Paczyriski (1981). 
Patterson (1998) favoured a modest increase of J on grounds 
of space density considerations and the ratio of systems be- 
low and above the period gap. Standard population models 
typically give a local CV space density of up to 10^* pc^"' 
(de Kool 1992, Kolb 1993, Politano 1996), with 99% of aU 
CVs below the period gap, 70% of these past period bounce. 
Patterson argues that the observed space density is at least 
a factor 20 smaller, and that there is little evidence for the 
predicted large population of post-period minimum CVs. 
A higher transfer rate after period bounce would remove 
systems from the population as the donor can lose all its 
mass in the age of the Galaxy, thus resolving both problems 
and the Pmin mismatch. However, postulating an as yet un- 
known J mechanism which conspires to produce almost the 
same value as Jgr at the transition from non-degenerate 
to degenerate stars does not seem very attractive. Rather, 
it seems more likely that the bulk of short-period systems 
are indeed unobservable - at least in the optical. Watson 
(1999) shows that the presence of a population of CVs with 
a space density of 10~* pc~^ that emits 10% of its accre- 
tion luminosity in X-rays cannot be excluded from ROSAT 
and ASCA data. It should be possible to place tight lim- 
its on such a population and its potential contribution to 
the Galactic Ridge emission from deep XMM and AXAF 
surveys. 

As a final note, we point out that the evolutionary ef- 
fect of an irradiation-induced stellar wind from the donor 
(considered semi-analytically by King & van Teeseling 1998) 
could help to resolve the Pmin mismatch and the period spike 
problem; detailed investigations are under way (Kolb et al. 
in preparation). 



5 SUMMARY 

We have investigated the problem that the standard expla- 
nation for the CV period minimum as a result of period 
bounce of systems driven by gravitational radiation predicts 
a period minimum that is too short, and an accumulation 
of systems close to Ptum, the "period spike", which is not 
observed. 

Using up-to-date stellar models by Baraffe et al. (1998) 
which successfully reproduce observed properties of single 
low-mass stars and brown dwarfs we confirm that Ptum is 
about 10 mins. shorter than the observed Pmin. 

We have presented synthesised period histograms for 
CVs below the period gap which, for the first time, are based 
on evolutionary calculations obtained with full stellar mod- 
els, and include CVs which form with brown dwarf donors. 
The period spike is always present. 

Although there are ways to explain why the spike is 
not observed — small number statistics for polars, unde- 
tectability for dwarf novae — we find no satisfactory reason 
why Pmin is louger than Ptum for both magnetic and non- 
magnetic CVs. 
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